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Abstract 

In this note, we present an algorithm that yields many new methods for constructing doubly stochastic and 
symmetric doubly stochastic matrices for the inverse eigenvalue problem. In addition, we introduce new 
open problems in this area that lay the ground for future work. 

Keywords: doubly stochastic matrices, inverse eigenvalue problem 
2000 MSC: 15A12, 15A18, 15A51 



1. Introduction 

An nx n matrix with real entries is said to be nonnegative if all of its entries are nonnegative. An nx n 
matrix A over a field F (F is either the real line R or the complex plane C) having each row and column 
sum equal to r, is said to be an r-generalized doubly stochastic matrix. The set of all n x n r-generalized 
doubly stochastic matrices with entries in F is denoted by fl r (n, F). A generalized doubly stochastic matrix 
is an element of f2(n, F) where 

n(n,F)= |J fl>,F). 

r£F 

Of special importance are the nonnegative elements in f2 1 (n, R) which are called the doubly stochastic 
matrices. The theory of doubly stochastic matrices is particularly endowed with a large collection of appli- 
cations in other area of mathematics and also in other disciplines (see for example (ll.l2.l^.l6lM.[l9ll25ll28l|). 
Let the set of all n x n doubly stochastic matrices be denoted by A„ and the set of all n x n symmetric 
elements in A„ will be denoted by A* . In addition, let M(n, F) be the algebra of all n x n matrices with 
entries in F and GL(n, F) be the general linear group over the field F. Finally, I is defined as the imaginary 
unit and for any matrix (or vector) A, its transpose will be denoted by A T . 

The Perron-Frobenius theorem states that if A is a nonnegative matrix, then it has a real eigenvalue 
r (that is the Perron-Frobenius root) which is greater than or equal to the modulus of each of the other 
eigenvalues. Also, A has an eigenvector x corresponding to r such that each of its entries are nonnegative. 
Furthermore, if A is irreducible then r is positive and the entries of x are also positive (see [1 IE EES HI). 
For doubly stochastic matrices, r — 1 and x = e n = {1,1, X) T ■ 

An intriguing object of study in the area of matrix theory and mathematical physics is that of the spectral 
properties and inverse eigenvalue problems for special kinds of matrices. The nonnegative inverse eigenvalue 
problem (NIEP) is the problem of finding necessary and sufficient conditions for an rt-tuplcs (Ai, . . . ,X n ) 
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(possibly complex) to be the spectrum of an n x n nonnegative matrix A. Although this inverse eigenvalue 
problem has attracted a considerable amount of interest, for n > 3 it is still unsolved except in restricted 
cases. Generally, we have two cases that result in three problems. 



• When A = (Aj., X n ) is complex, little is known. The case n < 3 was completely solved in [17( and 
the solution of the 4x4 trace-zero nonnegative inverse eigenvalue problem was given in (23[. 

• When the n-tuples (Ai, . . . , A n ) are real, then we have the following two problems: 

1) The real nonnegative inverse eigenvalue problem (RNIEP) asks which sets of n real numbers occur 
as the spectrum of an n x n nonnegative matrix A. 

2) The symmetric nonnegative inverse eigenvalue problem (SNIEP) asks which sets of n real numbers 
can occur as the spectrum of an n x n symmetric nonnegative matrix A. 

Each problem remains open. Many partial results for the three problems are known see the recent book [h| 
and the references therein, for the collection of all known results concerning these problems. Various people 
raised the question whether the (RNIEP) and (SNIEP) are generally equivalent. In the low dimension n < 4, 
the two problems are actually equivalent (see Q). However, for n > 4, the paper [l4| showed that the two 
problems are generally different. In addition, the paper [§| gives a construction for data A = (Ai,...,As) 
which is a solution for the (RNIEP) and there is no symmetric nonnegative 5x5 matrix with spectrum A. 

Another object of study in this area that has a big interest is the inverse eigenvalue problem for nonneg- 
ative matrices with extra properties. For example, we can consider the same problems for doubly stochastic 
matrices. So that we have the following problems. 

Problem 1.1. The doubly stochastic inverse eigenvalue problem denoted by (DIEP), is the problem of de- 
termining the necessary and sufficient conditions for a complex n-tuples to be the spectrum of an n x n 
doubly stochastic matrix. Equivalently, this problem can also be characterized as the problem of finding the 
region 0„ of C" such that any point in 0„ is the spectrum of an n x n doubly stochastic matrix. 

Now, when the n-tuples (Ai, . . . , A„) are all real, then we have the following two problems: 

Problem 1.2. The real doubly stochastic inverse eigenvalue problem (RDIEP) asks which sets of n real 
numbers occur as the spectrum of an n x n doubly stochastic matrix. This problem is also equivalent to 
the problem of finding the region 0£ of R n such that any point in Q r n is the spectrum of an n x n doubly 
stochastic matrix. 



Problem 1.3. The symmetric doubly stochastic inverse eigenvalue problem (SDIEP) asks which sets of n 
real numbers occur as the spectrum of an n x n symmetric doubly stochastic matrix. Equivalently, this 
problem can also be characterized as the problem of finding the region 0* o/R™ such that any point in 0* 
is the spectrum of an n x n symmetric doubly stochastic matrix. 

Example 1. The point a = (1,-1/2 + 7^3/2, -1/2 - is in 03 as a is the spectrum of A = 

( 1 \ 

On the other hand, the point A = (1,1/2,1/4) is in 0g as A is the spectrum of B = 



1 



V 



( 7/12 1/6 1/4 \_ 
1/12 2/3 1/4 
1/3 1/6 1/2 



V 



Moreover, A is in 0| as A is also the spectrum of C 



/ 13/24 7/24 1/6 \ 
7/24 13/24 1/6 
1/6 1/6 2/3 



V 



Although doubly stochastic matrices have been studied extensively, the (DIEP) and (RDIEP) have been 
considered in 22, 2(| where all the results obtained are partial. The (SDIEP) was studied in [ll| EM 2il 21, 
24 1, and earlier work can be found in [l3, 22, 27 1 and all the results obtained are also partial. For general 
n, all three problems remain open. In addition, the first and last of these problems have been completely 
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solved for n — 3 in 22]. Now in a similar manner to the (NIEP), one could raise the question whether the 
(RDIEP) and (SDIEP) are generally equivalent. For n = 3, we shall prove below that the two problems are 
actually equivalent. While for n > 4, it remains a very interesting open problem. In this paper, we are only 
able to obtain partial solutions concerning all 3 problems. 

Recall that X 6 f2 r (n, F) if and only if Xe n — re n and e^X = re^ if and only if XJ n = J n X = rJ n , 
where J„ is the n x n matrix with each of its entries is equal to —. Hence for any X 6 fl r (n, F), e n is an 
eigenvector for both X and X T corresponding to the eigenvalue r. In addition, for any V in GL(n, F) such 
that the first row of V is equal to a multiple of e„, and each of the last n — 1 row sums is zero, then clearly, 
its inverse V -1 has its first column equal to a multiple of e„ and each of the last n — 1 column sums of V^ 1 
is zero. Such V is said to have pattern S. Note that any n x n matrix V which is orthogonal and has its 
first column -^e n has pattern S. The main results in this paper rely on the following observation for which 

the proof can be found in [2(| ■ 

/ r \ 

Observation 1.4. Let V e GL(n, F) has pattern S and X e M(n - 1, F), then V^ 1 X \ V = A e 

f2 r (n, F). Conversely, for any A G tt r (n,F) and any V that has pattern S, there exists X G M(n — 1,F) 
/ r 

such that VAV- 1 = X 

In addition, an obvious necessary conditions for all three above problems are 



EiLi — f° r a U natural number k which just means that the trace of the nonnegative matrix A k 
is nonnegative. 

• < | Ai| < 1 for i = 2, n; as the Perron-Frobenius theorem insures. 

• If Xi is a complex eigenvalue of a doubly stochastic matrix D with nonzero imaginary part, then its 
conjugate is also an eigenvalue of D. 

Finally, we end this section with the following notation which is used throughout this paper. Let A be the 
n x n diagonal matrix with diagonal entries 1,A2,...,A„ with 1 > A2 > ... > A n > —1 and trace(A) = 
1 + A 2 + ... + A„ > 0. 



2. Construction of Matrices for the SDIEP 

In this section, we describe a way of obtaining partial solutions for the (SDIEP). As mentioned earlier, 
this problem is completely solved only for n = 3. Indeed the following theorem is proved in [22j. 

Theorem 2.1. 122] There exists a symmetric 3x3 doubly stochastic matrix with spectrum 1, A, \x if and 
only if -1 < A < 1, -1 < fj, < 1, A + 3^ + 2 > and 3A + n + 2 > 0. 

For n — 4, a conjecture is given in [TBI ] and for n > 4, only partial solutions are known. Now recall that 
is a convex polytope of dimension \n{n — 1), whose vertices were determined in 16( ( see also 



where it has been proved that if A is a vertex of A*, then A — l(P + P T ) for some permutation matrix 
P, although not every i(P + P T ) is a vertex. In addition, e„ is always an eigenvector of any n x n doubly 
stochastic matrix A so that when A is symmetric, then it is written as A = VqAVJ' for some Vo which is 
orthogonal and has pattern S. Therefore for a fixed such Vo one can ask what relation the {Xi} should 
satisfy for VoAV^f to be symmetric doubly stochastic. Then the region obtained in this way is a convex 
region of R" and the largest possible subregion of 9* one could obtain in this fashion, is attained when the 
columns of the matrix Vq are the common set of eigenvectors of a maximum number of vertices of A* that 
mutually commute (note that this is the case in 0, [53]). However this only gives partial solutions since 
the classification of all such Vq seems to be a difficult problem. Note that classifying all such Vq is equivalent 
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to the problem of finding the collection T> of all vectors in R n that can serve as the set of eigenvectors of 
an n x n symmetric doubly stochastic matrix. Since e n is always an eigenvector which is orthogonal to all 
other eigenvectors, then the collection S\{e„} is contained in the hyperplane of R" which is orthogonal to 
e n . Equivalently, if e„, a: 2 , x n are the orthonormal eigenvectors of a symmetric doubly stochastic matrix 
A corresponding to the eigenvalues 1, A2, A„ respectively, then A = ^e n e^ + X2X2X2 + ■■■ + X n x n Xn — 
J n + X2X2X2 + ■■■ + \ n x n Xn- Therefore, if x%, ...,x n are n — 1 vectors in R" such that the sum of the 
components in each of them is zero, then giving 1 > A2 > ... > X n > — 1, one can look at the conditions 
for which the matrix J n + \2X2X2 + ■■■ + X n x n x^ is nonnegative. Next, we state the earliest result for the 
(SDIEP) which is found in 



Theorem 2.2. ]Mj If 1 > A 2 > ... > A„ > -1 and 

11. 1 



-A 2 + 



-A3 + ...+ 



1 



-A„ > 



n n(n-l) (n-l)(n-2) "' (2)(1) 
then there is a symmetric doubly-stochastic matrix D such that D has eigenvalues l,A2,-..,A r , 



(1) 



G. Soules [27[ generalized the above result by considering the following orthogonal pattern S matrix 
(also known in the literature as a Soules matrix): 
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v / («- 


-l)(n-2) 
1 


\/n 
1 


y/n(n-l) 


v / («- 


-l)(n-2) 
1 


spn 
1 


y/n{n-l) 




-l)(n-2) 
1 


\fn 


y/n(n-l) 




-l)(n-2) 


1 


1 




1 


\fn 
1 


y/n(n-l) 




-l)(n-2) 
n-2 


\/n 
1 


y/n(n-l) 

71— I 




-l)(n-2) 




sfn 


^/n(n-l) 





1 
1 

-2 
v'O 





-1 

V2 






V 

It is easy to check (see (27{) that the symmetric matrix A - 
while the ith diagonal entry of A is given by the convex sum 



V n AVn has nonnegative off diagonal entries 



1 



E 



1 



)A„-fe+2 + ( — : — )A n _i+2, 



for i = 1, n. Moreover the an are increasing so the smallest one is an. So that if T is the convex region 
defined as the set of all n-tuples (Ai, A„) satisfying: 1>A2>...>A„>— 1 and 1 + A 2 + ... + A n > 
with an > for i = 1, n, then each point in T is a solution for the (SDIEP). Next, Soules obtains a larger 
region than T by constructing another matrix Vp from the matrix V n above as follows: 
Given an s-long strictly increasing sequence /3 with values in {1, ...,n} where < s < n, let /3 be the i-long 
sequence complementary to /3 (t = n — s). Define the s-vector u — (-^, ~^) T G R- s an d the t-vector 

u = ^) T € R n_;i . Now let w = (yf^)u and u) = -(^^)u. Finally, define B u as the s x (s - 1) 

matrix obtained from V s by deleting the first column and also B a is defined similarly. Next, for any /3 define: 

u w B u Oi 
Vp = I u w 2 Bn 



where Oi and 2 are respectively the s x (t— 1) and t X (s— 1) zero matrices. Then it is a routine computation 
to check (see 27]) that Vp is orthogonal and the symmetric matrix A(/3) = VpAVg has nonnegative off 
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diagonal entries. Now, if we let Vp be the matrix obtained from Vp by permuting its columns according to 
a permutation a that leaves columns 1 and 2 fixed and retains the relative order of the columns in B u and 
B u i.e. 

<jj = j for j = 1,2 and ctj < <Xj+i for j ^ s + 1, (2) 

then again the symmetric matrix A(cr, 0) = V» k(Vg) T would have nonnegative off diagonal entries. More- 
over, if a is an (s — l)-long strictly increasing sequence with values in {3, ...,n} defined by: a(j — 2) = <jj 
for j — 3, s + 1, and a is the (t — l)-long sequence complementary to a defined by: a(j — s — 1) = <jj 
for j = s + 2, ...,n, then each diagonal entry of A(a,f3) is either equal to E a or S a where E Q denote the 
following convex sum: 

£ = - + - ~ - A 2 + V Xa(s ~ k) 
a n ns ^-^ (k + l)fc 

and S a is obtained from E Q by replacing a with a. Finally, let n = 2m + 2 for n even and n = 2m + 1 for 
n odd. Now taking a to be the m-long sequence {3, 5, 7, n — 2, n}, then it is easy to see that E a can be 
obtained from S Q by replacing A„_2fc+2 with A n _2fc+i for k = 1, ...,m. So that S a is the smallest and then 
we have the following corollary: 

Corollary 2.3. 0/ // 1 > A 2 > ... > A„ > -1 and 

1 n-m-1 y> A n - 2 fe+2 r v 

fc— l 

holds, where n = 2m + 2 /or n even and n = 2m + 1 for n odd , then there exists an n x n symmetric doubly 
stochastic matrix D such that D has eigenvalues 1, A2, A n . 

2.1. Remarks on Soules's results. 

We can improve on the above corollary by choosing a 'better' permutation than a above. By that we 
mean a permutation that allows us to obtain a larger region of 0^, although this happens on the account 
of not having a unified condition for all cases of n. Indeed, we consider the following 4 cases: 

• n = 4fc + 1 : 

Taking a — {3, 4, 7, 8, n — 2, n — 1} and a = {5, 6, 9, 10, n — 4, n — 3, n} then we obtain the region 

(4) defined by: 

i-m-l \ , A 3 A 4 . A7 , a ti-2 , A n-1 - 

i ( m + 1 ) m ( m + 1 ) (m 1 ) m (m 2 ) ( m 1) " 6 + 2 ^ U , ^ 

1=T A 2 + ( m _l )m + (m-2)(m-l) + (m-3)(m-2) + " H 12 1 6 ' 2~~ — 

• n = Ak + 3 : 

let a = {3, 4, 7, 8, n — 4, n — 3, n} and then a = {5, 6, 9, 10, n — 2, n — 1} and we obtain the region 

(5) defined by: 

-1 \ I A 3 A4 _^ A 7 a ti-4 , A n-3 , An_ - „ 

1 ) m(m + l) "T" (m-l)m "T" (m-2)(m-l) "T" ■ ■ T 1 2 T 6 ~<~ 2 - U , K \ 

A„_ 9 , X„_i ^ „ W 



1_ I Tl-TTl X I A 5 _^ Afi , A9 . . A ?1 -2 , a ti-1 ^ n 

n T rem ^ 2 T (m-l)m T (m-2)(m-l) (m-3)(m-2) 6 T 2 — U 

n = 4fc + 2 : 

Take a — {3, 6, 7, 10, 11, n — 4, n — 3, n) and a = {4, 5, 8, 9, n — 2, n — 1} and the region (6) 
obtained is : 

{J_ , n-m-l \ , A 3 A 6 Aj^ A n-4 , A „_3 A „ 

71 ti(tti + 1) m(Tn + l) (m 1)tti 2)(?ji 1) "' 12 T 6 + 2 ^ U ff;) 

1 , n-m-1 , A 4 A 5 Ag A re-2 , A ..-l > n W 

TT + re(m + l) A2 + m(m + l) + (m-l)m + (m-2)(m-l) + '■■ 6 1 2^ ^ U 



5 



n = Ak + 4 : 

Here a = {3, 6, 7, 10, 11, n — 2, n — 1} and a = {4, 5, 8, 9, n — 4, n — 3, n} and the region obtained 
(7) is given by : 



7T + A2 + m(m + l) + I^lJ 



x 4 



(m-l)r, 



(m-2)(m-l) ^ 

+ (771-2)^771-1) + ' 



> 



> 



(7) 



Inequality (3) can be rewritten as 



n(rn + 1) 



A 2 + 



As 



m(m + 1) (m — l)m (m — 2)(m — 1) 



for n even 



n n{m + 1) 



A 2 + 



m(rn + 1) (n 



■ 2)(r. 



> 



> 



(8) 



0) 



for n odd. 

Since the {Aj} are in the decreasing order, then the left-hand side of (8) is less than or equal each of the 



left-hand sides of both (6) and (7). Also, note that ^- 



+ 



therefore the left-hand side 



n(m+l) m(m+l) 

of (9) is less than or equal each of the left-hand sides of both (4) and (5). So that the regions (4), (5), (6) 
and (7) are bigger. Here is an example that shows that the region (8) is a proper subset of the region (6). 
Consider (1, 1, 1, —1/2, —1/2, —1) e R 6 which clearly satisfies (6) and not (8) for n — 6. 

2.2. New Methods for Constructing Symmetric Doubly Stochastic Matrices 

Here we describe an algorithm for practical purposes that yields many new sufficient conditions for the 
symmetric doubly stochastic inverse eigenvalue problem. In fact, this can thought as a generalization of 
the construction described above. In addition, this algorithm can also be used with minor changes to find 
solutions for the (DIEP) in the complex case (see below). Before exploring this algorithm, it helps to think 
of the following simple observation which is the building block of the algorithm. 

Observation 2.4. Let x be an eigenvector of an n x n doubly stochastic matrix A corresponding to the 



eigenvalue A. // P denote the p x 1 zero vector, then 







and 



x 

0„ 



are respectively eigenvectors of the 



(n+p) x (n + p) doubly stochastic matrices I p © A and A © I p corresponding to X. 

Algorithm 1. Setpl: For low dimension k = 2,3,4,5,6, ... we consider all the vertices o/Aj*. For a chosen 
one of these vertices, we find its eigenvectors (using Maple for example) as orthonormal columns vectors 

Setp2: From the computed eigenvectors in Stepl, we form the k x k orthogonal pattern S matrix = 
(-jj=ek\x2\..-\xk) obtained from these eigenvectors, and if we let Afe be the kx k diagonal matrix with diagonal 
entries 1, A2, with 1 > A2 > ... > > —1, and 1 + A2 + ... + A^ > 0, then we check the conditions such 
that XfeAfcXj > for which by this construction, we always have solutions. 

SetpS: Next, for n > k we construct an n x n matrix W n by taking the first n — k + 1 columns of the Soules 

( X > \ 



matrix V n and the last (k — 1) columns 







( x k \ 




as the n columns of W n . Then find the region 



\ / V j 

for which W n AW^ is nonnegative. 

Setp4: Here we start the improvement process as follows. We construct another matrix Wp which will be 
the analogue of the matrix Vp (defined above) in the following manner: For k < s < n, let u be the s-vector 

(-7=, -7=), u be the (n — s)-vector (-7=, -7=), 



ir = ( ; ^j ± )u and w = — (y ^r^)^ be defined as in the 
previous section. Now define W u be the (s — 1) x s matrix obtained from the s x s matrix W s by deleting the 



first column and B a be the (n — s — 1) x (n — s) matrix obtained from the Soules matrix V n . 
the first column. Then we distinguish between the two cases: 



by deleting 
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Case 1: For n = k + 1, we iafce W n = Wk+i and we then find the region for which Wk+iXW k+1 is 



nonnegative. 
• Case 2: For n > k + 1, we let 



u w Oi W u 
Wp = | u w B a 2 



Oi and O2 are zero matrices of suitable orders. Then for any permutation a of the columns of Vp 
which preserves the relative order of columns in W u and P>u, the conditions for which WpAWj is 
nonnegative give a convex region T a C 0* whose all points are solutions for the (SDIEP). Finally, it 
is worth mentioning here that the union of all T a for all such permutation a, gives a union of convex 
subsets in Q s n whose again all points are solutions for the (SDIEP). 

Setp5 Repeat the same process for a different chosen vertex of A| . 

Remark 2.5. It should be noted that the above algorithm can also be used by starting with any k x k doubly 
stochastic matrix (i.e. any point of A s k ) that is not necessarily a vertex of A? which is the version of 
Algorithml that we will use for the case of (RDEIP) (see Section 3). However, the advantage of taking a 
vertex lies in obtaining a larger subset of Q s k which in turn by exploiting StepS and Step4 of Algorithml gives 
a larger region of Q s n . In addition, the above algorithm can also be used not just for low dimension, but for 
any dimension k for which there exists a k X k permutation matrix (or any k x k doubly stochastic matrix) 
P such that the eigenvectors of P can be computed. 



Remark 2.6. With Algorithm 1, the Soules matrix V n can be obtained by choosing the vertex P — 
of A| in Stepl since (t^ 5 "7f) T an< ^ ("71' — "7f) T are ^ e e ^9 envec ^ ors of P. 



Example 2. We apply Algorithm 1 on the matrix X = \ 



( 1 
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which is the orthogonal 



pattern S matrix that diagonalizes all the 4x4 zero-trace symmetric doubly stochastic matrices (see 
Note that the columns of X are the orthonormal eigenvectors of the following vertices of A|: 
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A simple matrix multiplication shows that XA.^X T = 
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Since the eigenvalues of A4 are in the decreasing order and trace(A4) > 0, then the only condition for which 
is nonnegative is simply given by : 



Ai - A 2 - A 3 + A 4 > 0. 



(10) 
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-1) 







case, 


is g 


iven by: 


1 


1 


1 1 


V20 


2 


2 2 


1 


-1 


-1 1 


V20 


2 


2 2 


1 


-1 


1 -1 


V20 


2 


2 2 


1 


1 


-1 -1 




2 


2 2 


-4 




































A = 



Then a simple matrix multiplication shows the matrices A — W n AW^ and B = V n AV^ differ only by 
the 4x4 principal submatrices formed from the first 4 rows and the first 4 columns of A and B i.e. 

x c \ , „ / y c 

C T D ) \C T D 

entries of the symmetric matrix A = W n AW^ are nonnegative except for the diagonal entries and for 

Ol2 = ^21 = a.34 = »43 = 

1 A2 A3 , A 4 , A„_3 A„_2 

n 



where X and Y are 4x4 matrices. So that by [27j all the 

tT 



An-l , A, 



n(rc-l) (n-l)(rc-2) (n-2)(n-3) 
In addition, the first 4 diagonal entries are equal to: an = 0,22 

_j_ A2 A3 r A 4 

n 



20 4 

a 33 = a 44 = 

, A„_3 A„_2 



A n - 



n(rc-l) (n-l)(rc-2) (n-2)(ra-3) 20 4 

and the remaining diagonal entries 055, a nn are increasing so the smallest one is 



An 

T' 



1 



A, 



A, 



A,, 



a 55 — — 

n 



We can rewrite 



_ 1 

0-55 — — 

n 



n(n — 1) (n — l)(n — 2) 

A2 A3 A„_3 

20 



30 

An-3 



16 A 



A n - 



A n - 



n(n — 1) (n — l)(n — 2) 

As 1 > A2 > ... > A„ > — 1, then clearly an < 055 and therefore the conditions for which A = W^[AW n is 
nonnegative are: 

'an > 

flu >0 (11) 



Hence we have the following theorem. 

Theorem 2.7. Let 1 > A 2 > ... > A„ > -1 arcd 1 + A 2 + ... + A„ > 0. // 

( 1_ 1 A2 1 A 3 , Xa , , A„_ 3 , A„_ 



i(n-l) 
A 2 



(n-l)(n-2) 
A3 



(n-2)(n-3) 
A 4 



20 



4 

An-2 



An-l 
4 

A„_i 



^ >o 

Y >0 



n T ra(n-l) ' (n-l)(n-2) r (ra-2)(n-3) T *" ' 20 4 4 

i/ien i/iere is an n x n symmetric doubly stochastic matrix D such that D has eigenvalues 1, A2, A ri 



u w Oi W u 

To improve upon the results of the above theorem, we consider the matrix Wp = I u id B a O2 

Now let a be the permutation of the columns of Wp (which of course preserves the relative order of columns 
in W u and B a ) given by: 

• For n = 2m + 2 even and s — m + 1 and a = {4, 6, 8, n — 6, n — 3, n — 1, n} then the new improved 
conditions become: 



^1 X, J- A 4 1 A 6 j ^8 1 1 A "-6 1 A n-3 , A n-1 , A, ^ » 

FTJ 2 "T" m(m + l) (m-l)m "T" (m-2)(m-l) T ■■■ "I" 2 4 4 ~T~ 4 — U 

X„ J- A 4 I A 6 j ^8 L I A "-6 _ A n-3 _ A n-1 i A" > n 

. 1( ^2~T ^(^ + 1) ' (tti — l)m ' (??i-2)(??i-l) i ■■■ I > ( i i i 



20 1 


■) ■ .) 
4 




X n _ 6 X. 


1 - 3 


A „-i 


20 


1 


4 


n — 6, n - 


-3. 


n— 1, 


*n-6 , \ 


-1-3 


a „-i 


20 1 


4 


T 4 


A„_ 6 A. 


a - 3 


A n-1 



(12) 



• wnue tor n = zm + i odd and s = m + 1, and a = \6, o, / 
conditions become: 

J_ , n-m-1 , , A 3 I A 5 j *7 I I ±L 1_ : _j 1 - 1 _1_ An, > n 

n "I" „( m + l) A 2 -T m(m + 1) -T ( m _i) m -I- ( m _2)(m-l) T 20 4 T 4 T 4 ^ u ,, „N 

1 , n-m-l , A3 A 5 A 7 A „- 6 A„_ 3 A „_l A„ > „ ' ' 

l> "T n(m + l) A 2 "I" m ( m + i) T ( m _l) m T ( m _2)(m-l) ■ ■ ■ 20 4 4 ~i 4 _ u 

3. The Real Inverse Eigenvalue Problem For Doubly Stochastic Matrices 

First as mentioned earlier, this problem has been considered in [26[ where the following theorem has 
been obtained. 

Theorem 3.1. \2a l Let a = {1, A2, A„} be a set of real numbers such that 

1 > A 2 > ... > X r > > A r+ i > ... > A„. 
If 1 > A2 + nmaa;{|A2|, |A n |}, then there exists an n x n doubly stochastic matrix with spectrum a. 

In addition, the case n = 3 has been used in a proof of a theorem concerning the (SDIEP) in 22] in which 
the following result holds. 

/ P q l-p-q \ 

Theorem 3.2. Let X = I r s 1 — r — s \ be a doubly stochastic matrix with 

\ 1 — P — r 1 — q — s p + q + r + s — 1 J 
real eigenvalues 1, A, fi. Then —1 < A < 1, —1 < < 1, A + 3/i + 2 > and 3A + fx + 2 > 0. 

Proof. The first two inequalities are obtained from the Perron-Frobenius theorem. For the other two 
inequalities, let x = A + 3/i + 2 and y = 3A + fi + 2. We prove that x and j/ are nonnegative by showing 
that their sum x + y and their product xy are nonnegative. For, x + y = 4A + 4// + 4 = Atrace(X) > 0. 
Now z?/ = 3(A + ^ + l) 2 + 2(A + /i + 1) + 4A/i - 1 = 3 [trace (X)} 2 + 2trace(X) + Adeterminant{X) - 1 = 
3(5 — r) 2 + 12(p + s)(p + q + r + s — 1) + 12ps. As the entries of X are nonnegative, therefore xy > 0, and 
the proof is complete. ■ 

Combining the above theorem with Theorem 12. 1[ we have the following conclusion. 

Corollary 3.3. The (RDIEP) and the (SDIEP) are equivalent for the case n = 3. 

Now returning to the procedures described at the beginning of Section2 and starting out with a particular 
n x n pattern S matrix V which is not orthogonal, we obtain the following theorem that solves the (RDIEP) 
for some restricted cases. 

Theorem 3.4. If 1 > A 2 > ... > A„ > -1 and 

l-(n-l)A 2 + A 3 + ... + A„>0 ( . 

l + (n-l)A„>0 1 ' 

then there is an n x n nonsymmetric doubly stochastic matrix D with spectrum 1, A2, A„. 

9 



Proof. Clearly the matrix 



/ 1 



V = 



V 1 



1 i 

-1 

-1 











-1 



1 







-1 





1 \ 










-1 / 



has pattern S and its inverse V 1 is given by: 



/ 



(n-l) 



1 

n 
1_ 

n 

(n-1) 



(n-1) 



V h 



1 

n 

(n-1) 



Now the entries of the matrix A = (a^ ) = VAV 1 satisfy the following relations: 



an 

an - 

an : 



L)A< + Xi+i + ... + A„) 



(15) 



i(tmce(A)) 

i(l + (n-l)Ai) fori = 2,...,n 
i(H-A 2 + ... + A i _i-(n-l)A i 
|(1 - Aj) for j > 1 and j ^ i 

Note that an and djj = ^(1 — Aj) for j > 1 and j ^ i are nonnegative since the diagonal entries of A are 
in the decreasing order. In addition, for i ^ 1, the entries an are increasing so the smallest one is 021 and 
an are decreasing so the smallest is a nn . Therefore the matrix A is nonnegative if and only if ai\ > and 
a nn > 0. Finally A is doubly stochastic since V has the pattern S, and the proof is complete. ■ 

It is should be noted here that the conditions of the above theorem differ than those of Theorem 13.11 as 
the point (1, 1/2, 1/4) clearly satisfies the conditions of the above theorem but obviously does not satisfy 
the conditions of Theorem 13. II for the case n = 3. In addition, the conditions of the above theorem are also 
sufficient for the existence of an n x n symmetric doubly stochastic matrix D with the spectrum (1, A2, A„). 
To see this, it suffices to look at the following result which is Corollary 7 in [l2^ . 

Theorem 3.5. flB] If X2, A n G [— l/(n — 1), 1], then there exists an n x n symmetric doubly stochastic 
matrix D with spectrum (1, A2, A„). 

Here it is worth mentioning that at this stage and with extensive numerical computations, we are not 
able to find a list of n (> 4) real numbers which shows the two problems (RDIEP) and (SDIEP) are different. 
In conclusion, the question whether (RDIEP) and (SDIEP) are equivalent or not, for n > 4 remains an open 
problem. 

We conclude this section by using the version of Algorithm 1 mentioned in Remark 2.5. More explicitly, 
we choose a nonsymmetric doubly stochastic matrix with real eigenvalues such as the 3x3 matrix B 
of Example 1. Using Maple, the eigenvectors of B are given by (1,1, 1) T , (1,— 2,1) T and (1,1,— 2) T . 

1 1 1 

Now the corresponding pattern S matrix is given by P = | 1 —2 1 | . Then its inverse is given by 

11-2 



10 



. Then the conditions 



1/3 1/3 1/3 
P- 1 = | 1/3 -1/3 
1/3 -1/3 

for which PA 3 P _1 is nonnegative is given by Theorem 3.4 for the case n — 3 with no surprise as P is a 
multiple of V^ 1 which is used in the proof of Theorem 3.4. Following Step4 of Algorithm 1, we construct 
the matrix W n which is given by: 



1+A 2 + A 3 


1-A 2 


I-A3 


3 


3 


3 


1-2A 2 +A 3 


1+2A 2 


I-A3 


3 

I+A2-2A3 


1-A 2 


3 

1+2A 3 


3 


3 


3 



( J_ 

s/n 
1 

s/ri 
1 



1 

s/ri 



1 

s/ri 

1 

s/ri 
1 



sjn(n- 


-1) 




-l)(n- 
1 


-2) 


Jn(n- 
1 


-1) 




-l)(n- 
1 


-2) 


y/n(n- 


-1) 


\A«- 


-l)(n- 
1 


2) 


\Jn(n- 


-1) 


>/("- 


-l)(n- 


2) 


1 






1 




sjn(n- 


-1) 


v / («- 


-l)(n- 

71-2 


-2) 


sjn(n- 

n — 


-1) 
1 




-l)(n 



-2) 


y/n(n 


-1) 







1 

V^L2 


1 


1 


1 

s/12 


-2 


1 


1 

v^L2 


1 


-2 


-3 
VT2 




































A virtually identical description as in the symmetric case shows that A = W n KW n 1 



B 



where X and Y are 3x3 matrices. Thus a simple check shows that the 



Y C 
C T D 

conditions for which A = WnAW,^ 1 is nonnegative give the following theorem. 



VnAV- 1 = 



X 
C T 



C 
D 



and 



Theorem 3.6. Let 1 > A 2 > ... > A„ > -1, such that 1 + A 2 + ... + A„ > 0. If 

2A 



A 2 



+ ■ 



n(n-l) " r (n-l)(n-2) 
I _|_ A 2 A 



+ 



(n-2)(n-3) + - + (4)(3) " 
n(n-l) + (n-l)(n-2) + (n-2)(n-3) + - + WW) + ~3~ - U 



A„ 



£/ien (1, A 2 , A„) is i/ie spectrum of an nx n nonsymmetric doubly stochastic matrix. 

Next, we use the improvement process of Step4 to obtain new results. This can be illustrated by the following 
example for the case n = 6. 

Example 3. we apply the improvement process described in Step4 of Algorithm 1 to obtain the matrix 



Wp 



± 
± 

s/S 
V 76 



1 

't 

s/6 









1 


1 


\ 








-2 


1 










1 


-2 




1 


1 










-2 


1 










1 


-2 








/ 



Now an inspection shows that WpAeWp 1 



B C 
C D 



where 



B 



1 I A2_ I _A_5_ _|_ A^ 

6 ~ r .6 ^ 3 ~ 3 

J_ 1 A^ 2A5 I Ae 

6 ' 6 ,3 ^ 3 

J_ 1 A^ _j_ _A_5_ 2A6 

6 6 3 3 



A2 
.6 



A 



A 2 _|_ 2A5 _L _|_ A2 A6 £1 



_ Ae 

A 2 Ae_ 



I -^2 



1 A2_ 1 ^A 6 



1 


A 2 


1 


A 2 


1 


A 2 


6 




6 




6 




1 


A 6 2 


1 


A 6 2 


1 


A 6 2 


6 




6 




6 




1 




1 


A 6 2 


1 


A 6 2 


6 


6 


6 


6 


6 


6 
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and L> = h + ¥ ~ ^¥ + ¥ I + ¥ + ^¥ 7 + ¥ - ¥ • Note that for 1 > A 2 > ... > A 6 > -1, 




A-2 


A 3 


1 


| A 2 


A 4 






6 


A 6 






, 2:^ 


1 


1 A2 

A 6 

+ ¥ 


_£ 




A 3 
A3 


6 
1 


. 2X4 


6 


3 


6 




r 3 



we have i + ¥ + ^ ^ | + if + 3 ■ Thus the conditions for which the matrix WpK^Wg 1 is nonnegative 
result in the following conclusion. 

Theorem 3.7. Let 1 > A 2 > ... > A 6 > -1, suc/i t/iat 1 + A 2 + ... + A 6 . // 



1 1 A2 _i_ 2Ae 

6 ' 6 ' 3 
L 1 A2 2 A3 I A4 



> 



6 1 6 3 1 3 — 

j_ 1 A2_ 2A5 I Ae_ \ p. 

6^6 3 ~t~ 3 - u 



then (1, A 2 , Xq) is the spectrum of a 6 x 6 nonsymmetric doubly stochastic matrix 



4. Constructing doubly stochastic matrices with complex spectrum 



We start the section by again mentioning that only the case n = 3 is completely solved in [22j where the 
following theorem has been proved. 

Theorem 4.1. f2$] Let z be a complex number with nonzero imaginary part and z be its complex conjugate. 
Then (1, z, z) is the spectrum of a 3 x 3 doubly stochastic matrix if and only if z is in the convex hull of the 
three cubic roots of unity 

To take advantage of Algorithm 1 for this case, we recall the famous Birkoff 's theorem that states that 
A„ is a convex polytope of dimension (n — l) 2 where its vertices are the n x n permutation matrices. 
Next, we describe how to manipulate Algorithm 1 to obtain solutions for (DIEP) with complex spectrum 
but with minor changes namely in Stepl we start out with a nonsymmetric vertex of A^. In fact, we can 
begin Algorithm 1 with any k x k doubly stochastic matrix A with complex spectrum but in this case the 
orthonormalization process in Stepl should be dropped if the columns of A are not orthogonal. We illustrate 
this idea by considering the following example. 

1 \ 

Example 4. The vertex P = | 1 I of A3 has the following three eigenvectors -¥^3, %2 = 

10 0/ 

-^(1, w 2 , w) T and X3 = ;4g(l) w, w 2 ) T where w = —1/2+/V3/2 is the primitive cubic root of unity. Then the 

11 1 

3x3 complex pattern iS matrix obtained from these eigenvectors is given by X = | 1 w 2 w | and its 



1 



:,2 



11 1 \ / 1 

inverse is A" 1 = ( 1 w w 2 . If we let II3 = a + lb ) with trace(Il 3 ) = 1 + 2a > 

' 1 w 2 w ) \ a- lb 

and I a + Lb\ < 1, then a simple check shows that 

I 1 1 + 2a 1 - a - by/3 1 - a + byfl 

AI^A" 1 = - ( l-a + 6V3 l + 2a 1-a-by^ 
1 - a - bVi 1 - a + bV3 1 + 2a 

Therefore the conditions for which AI^A -1 is nonnegative are simply giving by: 



1 - a + bV3 > 
1 - a - b^/3 > 



(16) 



12 



The next step of our algorithm is to form the following n x n complex pattern S matrix 



7™ 
1 


\/™(n- 
1 


-1) 




-l)(n- 
1 


-2) 


7™ 
1 


y/n(n- 
1 


-1) 




-l)(n- 
1 


2) 


7™ 
1 


^/n(n- 


-1) 


\/(n- 


-l)(n- 
1 


2) 


y/n 


yjnin- 


-1) 




-l)(n- 


2) 


1 


1 






1 




V" 
1 


\/™(n- 


-1) 




-l)(n- 
n-2 


-2) 


y/n 

1 


n — 


-1) 
1 




-l)(n 

n 


-2) 



1 


1 


1 


\ 




V o 


V 3 




1 








7l2 




73 




1 


to 


™ 2 




7i2 








-3 
7l2 









































/ 



\ 7™ ^/n(n-l) 

Now if we let II„ be the diagonal matrix 1 © A2 © ... © A„_2 © a + 76 © a — 76, such that |a + 7o| < 1 and 
with 1 > A 2 > ... > A„_ 2 > — 1 and trace(TI n ) = 1 + A 2 + ... + A„_ 2 + 2a > 0. Then an inspection shows that 
A B \ 

T ^ \ where A is a 3 x 3 circulant matrix whose first row (an, ai 2 , 013) is given by: 



an = 


j_ 

n 


+ 


ai2 = 


1 

n 


+ 


ai3 = 


1 

n 


+ 



n(n — 


17 + 


(n-l)'(n 


-2) 


(n-2)(n- 


-3) 


A 2 




A3 




A4 




n(n — 




(n-l)(n 


-2) - 


r (n-2)(n- 


-3) 


A 2 




A3 




A4 





A„-2 


, 2a 




(4)(3) 


~T 3 




A„-2 


a 


b73 


(4)(3) 


3 


3 


A„-2 


- - + 


675 


(4)(3) 


3 ^ 


3 



and 73 is nonnegative and C has nonnegative off-diagonal entries and its diagonal entries are increasing so 
that the smallest is 



Cll 



1 

- + 



A 2 



A, 



A 4 



A n - 



3A 



ra-2 



n n(n-l) (n-l)(n-2) (n-2)(n-3) "' (5) (4) 4 
Thus the conditions for which WrJlnW' 1 is nonnegative gives the following theorem. 
Theorem 4.2. Let 1 > A 2 > ... > A„_ 2 > -1, such that 1 + A 2 + ... + A„_ 2 + 2a > and \a + Ib\ < 1. If 

1_ 1 A 2 A3 1 A 4 I, A„_2 1 2a 

71 
J_ 

n 
1_ 

n 
1 



^2 

n(n — 


1) 


+ (n- 


^3 

l)(n- 


2) + 


"\4 

(n-2)(n 


-3) 


A 2 






A 3 




A 4 




n(n — 


1) 


+ (»- 


l)(n- 


2) + 


(n-2)(n 


-3) 


A 2 






A 3 




A 4 




n(n — 


1) 


+ (»- 


l)(n- 


2) + 


(n-2)(n 


-3) 


A 2 




+ 7— 


A3 


+ 


A 4 





(4)(3) 

A„-2 

(4)(3) 

A„-2 



f >0 



3 

3A„ 



6yJ > 
6yJ > q 



> 



■ (5)(4) ^ 4 

i/ien (1, A2, A n _ 2 , a + ^o, a — lb) is the spectrum of an n x n doubly stochastic matrix 

Note that the above theorem deals with the case of at most two nonreal eigenvalues namely a + lb and 
a — lb in the list {1, A2, A„_ 2 , a + Ib, a — lb}. Now, we can use the improvement process of Step4 to obtain 
conditions on a list with more than just two complex eigenvalues. This can be illustrated by the following 
example again for the case n = 6. 

Example 5. we apply the improvement process described in Step4 of our algorithm to obtain the following 
matrix 



( A 
1 

7e 
1 

\/6 
76 

1 

76 

\ Ve 



1 

V6 
1 

76 
1 

1. 

76 
1 
76 

1 

V6 








1 

73 








1 

73 



75 73 

w w 2 

73 73 



1 


1 


w 
73 




W 

73 


73 
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If wc let a be the diagonal matrix 1 (B f (B a + lb (B a — lb (B c + Id (B c — Id, such that \a + Ib\ < 1, 
|c + Id\ < 1, — 1 < / < 1 and trace(a) = 1 + / + 2a + 2c > 0. Then an inspection shows that the matrix 

WpaWp 1 has the form WpaW^ 1 = ( q ^ I where 



B 



6 ^ 6 T 
1 i/_ c 
6 ' 6 3 
1 i/_c 
6^6 3 



2c 
3 



3 

rfV3 
3 



, / _ c _ dy^ 
'63 3 

6 ~ 6 ~ 3 
, / _ c , dV3 
6 3 3 



/ _ £ 

6 3 

/ _ £ 

6 3 



_ cK/3 
3 

_ 

3 

2c 



1 _ / 1 _ / 1 _ / 

6 6 6 6 6 6 

c=l i - i i - i i - 2 

' w I 6 6 6 6 6 6 

!_/ I _ 1 I _ Z 

6 6 6 6 6 6 



and 



J_ i X a_ fr-s/3 1_ i / a_ _i_ fry^3 

6 ' 6 3 3 6 ' 6 3 " r 3 

h/A Ii/|2a 1 , / _ a _ b\/3 

3 6 ' 6 ^ 3 6^6 3 3 

a. fo-s/3 1 I / a_ i b-s/3 J_ _i_ / _j_ 2a 

3 3 6 ' 6 3 ' 3 6 " r 6 ' 3 

WpaWp 1 is nonnegative result in the following conclusion. 



6 ~ 6 ~ 3 
_|_ X _ £ 
'6 3 



4 — # — • Now the conditions for which the matrix 



6^6 



Theorem 4.3. Let 1, f,a + lb, a — Ib,c + Id,c— Id be complex numbers such that \a + Ib\ < 1, \c + Id\ < 1, 
-1 < / < 1 and 1 + / + 2a + 2c > 0. // 



I _|_ 1 _|_ 2a > Q 

6 T 6 T 3 — u 

I + L + 2£ > 

6 ~ 6 ~ 3 — u 

1 i / _ o i fc^3 -> pi 

6^6 3 ^ 3 — 

1 i / _ a _ fc\Z3 -v. p, 

6^ 6 3 3 — u 

I _i_ L _ £ _|_ rf V3 > 

6 6 3 ' 3 — u 

1 i / _ c _ dV3 

6^6 3 3 



> 



then (1, /, a + 76, a — 76, c + Id, c — Id) is the spectrum of a 6 x 6 doubly stochastic matrix 

Finally, it is worth mentioning here that even applying only the first two steps of our algorithm has 
an interest of its own since it has the advantage of yielding new conditions for the (DIEP) for the chosen 
dimension k in Stepl. To illustrate this, we include the following example. 



/ 

Example 6. Consider the following 4x4 nonsymmetric permutation matrix p = 

V 

Maple for example, the complex eigenvectors of p are given by: x\ = (1, 1, 1, 1), x 2 
{—1, 1,-1, /), and x 4 



1 





(V 



\ 



Using 



-1, -1, 1), £3 = 



X = 



V 



\ 



_L 

2 
_ 1 
2 

I 

2 



(1,1,-1, —I). Then, from these eigenvectors, we form the complex pattern S matrix 
1 

2 

_L 

2 
_ 1 
2 

_ I 

2 



In addition, wc let IT4 be the diagonal matrix with diagonal entries as 1, c, 



a + Ib, a — lb, with trace(n 4 ) = 1 + c + 2a > 0, —1 <c< 1 and \a + Ib\ < 1. Next, we check the conditions 
for which X^X -1 is nonnegative as follows. A simple matrix multiplication shows that 



xn 4 x- 



( \ + 



1 

4 

1 _ c 

4 4 

1 _ c 

4 4 

- + - 

4^4 



4 2 
c _i_ £ 
4^2 



4 1 2 

£ a 

4 2 

£ 1 a 
4^2 

£ _ £ 

4 2 



I + £ _ 
4^4 

I _ £ _ 
4 4 

I _ £ + 
4 4 ~ 

I + £ + 
4 4 



/ 



Note that each diagonal entry of XU4X 1 is equal to ^ trace (II4) and therefore is nonnegative. Hence the 
conditions for which XU^X^ 1 is nonnegative, are given by: 



I _ £ _ £ > n 
4 4 2 — u 

I - £ 4- £ > 
4 4 ^ 2 — u 
I 4- £ - £ > 
4^4 2 — u 

14 



(17) 



Thus we have the following theorem. 

Theorem 4.4. Let 1, c, a + Ib, a — lb, be complex numbers with l + c + 2a> 0, — 1 < c < 1 and \a + Ib\ < 1. 
If 

1 - c - 26 > 
l-c+26>0 , 
1 + c - 2a > 

£/ien (1, c, a + 76, a — 76) is i/ie spectrum of a 4 x 4 doubly stochastic matrix D. 

Note that if we wish to have the realizing matrix D has zero trace, then it suffices to add the extra constraint 
trace(Q) = l + c + 2a = and then we obtain the following. 

Theorem 4.5. Let 1, c, a + Lb, a — lb, be complex numbers with l + c+2a= 0, —1 < c < 1 and \a + Ib\ < 1. 
If 

1 - c - 2b > 
1 - c + 26 > , 
-1 < a < 

then (1, c, a + 76, a — 76) is the spectrum of a Ax A doubly stochastic matrix D with zero trace. 
Conclusion 

We described an algorithm that yields many sufficient conditions for three inverse eigenvalue problems 
concerning doubly stochastic matrices. Although at this stage we did not offer complete solutions, we leave 
it for future work to check if it is possible to characterize all the necessary points in Stepl needed for this 
algorithm to offer complete solutions at least for low dimensions. However, besides that it offers many new 
partial results, the main importance of this algorithm lies in the fact that it can be used as a checking point 
in case of a conjecture concerning complete solutions is given for any of these three interesting problems. 
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